Optimization Of Neutron-Gamma Tools For Inelastic-Gamma Ray Logging

ABSTRACT

Systems, methods, and devices for inelastic gamma-ray logging are provided. In one embodiment, such a method includes emitting neutrons into a subterranean formation from a downhole tool to produce inelastic gamma-rays, detecting a portion of the inelastic gamma-rays that scatter back to the downhole tool to obtain an inelastic gamma-ray signal, and determining a property of the subterranean formation based at least in part on the inelastic gamma-ray signal. The inelastic gamma-ray signal may be substantially free of epithermal and thermal neutron capture background.

BACKGROUND

The present disclosure relates generally to oilfield logging and, moreparticularly, to oilfield logging involving inelastic gamma-raymeasurements.

Many properties of a subterranean formation may be determined usingdifferent oilfield logging techniques, which may involve one or moretools having a radioisotope source. For example, to locate gas in asubterranean formation, a conventional practice combines data obtainedfrom two tools. One of the tools is a “density” tool, which measures theelectron density of the formation, and the other of the tools is a“neutron porosity” tool, which generally measures the density ofhydrogen in the formation, known as the “hydrogen index (HI).” Based onmeasurements of formation density and hydrogen index, the porosity andpore fluid density of the formation may be determined. For a givenformation fluid density, or gas saturation, a combination of a decreasein the formation density and an increase in the hydrogen index indicatesan increase in the porosity of the formation. Meanwhile, for a givenformation porosity, a combination of a decrease in the formation densityand a decrease in hydrogen index indicates a decrease in the pore fluiddensity and hydrogen content. For pores filled with water and gas or oiland gas, the density and hydrogen index are an indication of the gassaturation (volume fraction of the pores occupied by gas). For poresfilled with gas only, the density and hydrogen index are an indicationof gas density (pressure).

The density and neutron porosity tools for measuring formation densityand hydrogen index may generally employ radioisotope sources to obtainformation density and hydrogen index measurements, respectively. Forexample, the density tool may use a source such as ¹³⁷Cs to emitgamma-rays into a formation. Based on a count of gamma-rays scattered bythe formation, the density tool may determine the electron density ofthe formation. Similarly, the neutron porosity tool may use a sourcesuch as ²⁴¹AmBe to emit neutrons into a formation. A count of neutronsscattered by the formation may yield a hydrogen index measurement. Suchradioisotope sources may be disadvantageous in oilfield tools, as thesources may be heavily regulated by law and they can be hazardous sincethey cannot be shut off.

In lieu of such radioisotope sources, an electronic neutron generatormay be used which will produce neutrons which, in turn, producegamma-rays. To do so, the electronic neutron generator may emit neutronsinto a formation, which may in turn produce gamma-rays via inelasticscattering and neutron capture events. A count of gamma-rays produced byinelastic scattering may generally yield a signal that is related toformation density, and a count of scattered neutrons may generally yielda neutron porosity signal that corresponds to the hydrogen index of theformation. Alternatively, a count of capture gamma-rays may generallyyield a neutron porosity signal that corresponds to the hydrogen indexof the formation. If it is not possible to separate the inelastic andcapture gamma-rays to produce nearly independent signals sensitive toformation density and hydrogen index, respectively, then the two signalsmay not be used together to enable a precise determination of porosityand gas saturation.

Neutron reactions that produce gamma-rays may be separated according tothe energy of the neutron. After a 14 MeV neutron has been emitted bythe source, it begins to lose energy by the processes of elastic andinelastic scattering. Inelastic scattering events are typically producedby neutrons in the energy range 1-14 MeV. After neutrons have decreasedin energy below approximately 1 MeV, they typically have insufficientenergy to inelastically scatter; however, they continue to lose energyby elastic scattering. The decrease in energy from 14 MeV to 1 MeVhappens very rapidly, in a time typically less than 1 microsecond.Inelastic scattering reactions therefore occur very quickly after theneutron leaves the source, typically in less than 1 microsecond. Fromapproximately 1 MeV down to thermal energy (approximately 0.025 eV),neutrons decrease in energy by elastic scattering over a time intervalthat ranges from 2 to several microseconds, depending on the amount ofhydrogen in the formation. During that slowing time, neutrons may becaptured and this may lead to the emission of one or more gamma-rays.These are so-called “epithermal” capture gamma-rays. Neutrons whichdecrease in energy completely to thermal energy continue to elasticallyscatter at that energy, often for many hundreds of microseconds untilthey are captured and this may lead to the emission of one or moregamma-rays. These are so-called “thermal” capture gamma-rays. Sinceneutrons are emitted from an electronic neutron source typically inbursts no shorter than 10 microseconds, it will be appreciated that theinelastic and epithermal capture gamma-rays are emitted substantiallywithin that 10 microsecond burst and therefore overlap in time. Thermalcapture gamma-rays, on the other hand, extend into the time intervalbetween bursts as well as during the burst. Since there is overlap ofcapture and inelastic gamma-ray events during the burst, simply summingover all counts during the burst may yield a signal that corresponds, atleast in part, to the hydrogen index of the formation, rather than toformation density. Various techniques to correct for the component ofthermal neutron capture gamma-rays have been disclosed, for example, inU.S. Pat. No. 5,374,823 to Odom. Odom suggests that epithermal capturegamma-rays might be corrected if the epithermal neutron lifetime wereknown but gives no guidance on how to do this and suggests that it isunimportant anyway. Trcka, in U.S. Pat. No. 7,365,308 mentions theproblem of epithermal capture gamma-rays but is silent on how to correctfor them. In a similar way, Wilson in U.S. Pat. No. 6,207,953 discussesan “inelastic” gamma-ray count rate but this is just the burst sumcorrected for thermal capture gamma-rays. No attempt is made to correctfor the epithermal gamma-ray contamination. However, unless both thermaland epithermal capture gamma-rays are eliminated from the burst sum, theresulting signal will generally correspond in large part to hydrogenindex rather than formation density, and is redundant with a capturegamma-ray or scattered neutron measurement of hydrogen index.

SUMMARY

Certain aspects commensurate in scope with the originally claimedembodiments are set forth below. It should be understood that theseaspects are presented merely to provide the reader with a brief summaryof certain forms the embodiments might take and that these aspects arenot intended to limit their scope. Indeed, the embodiments may encompassa variety of aspects that may not be set forth below.

The present disclosure relates to systems, methods, and devices forinelastic gamma-ray logging. In one embodiment, such a method includesemitting neutrons into a subterranean formation from a downhole tool toproduce inelastic gamma-rays, detecting a portion of the inelasticgamma-rays that scatter back to the downhole tool to obtain an inelasticgamma-ray signal, and determining a property of the subterraneanformation based at least in part on the inelastic gamma-ray signal. Theinelastic gamma-ray signal may be substantially free of epithermal andthermal neutron capture background.

BRIEF DESCRIPTION OF THE DRAWINGS

Advantages of the presently disclosed subject matter may become apparentupon reading the following detailed description and upon reference tothe drawings in which:

FIG. 1 illustrates a system for inelastic scattering well logging, inaccordance with an embodiment;

FIG. 2 illustrates the total gamma-ray detector counts as a function oftime in accordance with an embodiment;

FIG. 3 illustrates a manner of operating the system of FIG. 1, inaccordance with an embodiment;

FIG. 4 is a plot illustrating thermal and epithermal neutron absorptioncapabilities of a lithium-containing shielding material, in accordancewith an embodiment;

FIG. 5 is a flowchart describing an embodiment of a method fordetermining porosity and gas saturation using the system of FIG. 1, inaccordance with an embodiment;

FIG. 6 illustrates the total gamma-ray detector counts as a function oftime in an embodiment using materials with a low resonance integral;

FIG. 7 is a flowchart describing an embodiment of a method for removingan epithermal neutron capture background from a gamma-ray count, inaccordance with an embodiment;

FIG. 8 is a flowchart describing another embodiment of a method forremoving an epithermal neutron capture background from a gamma-raycount, in accordance with an embodiment;

FIG. 9 is a plot illustrating an exemplary count of gamma-rays detectedduring and immediately following a neutron burst using the system ofFIG. 1, in accordance with an embodiment;

FIG. 10 illustrates the total gamma-ray detector counts as a function oftime in an embodiment using materials with a low resonance integral anda method of subtracting capture background;

FIG. 11 is a flowchart describing another embodiment of a method forreducing an epithermal neutron capture background from a gamma-raycount, in accordance with an embodiment;

FIG. 12 is a flowchart describing another embodiment of a method forreducing an epithermal neutron capture background from a gamma-raycount, in accordance with an embodiment;

FIG. 13 is a plot simulating the determination of porosity and gassaturation of a formation using the system of FIG. 1, in accordance withan embodiment; and

FIG. 14 is a flowchart with an alternate approach for determiningporosity and gas saturation without using a signal from a neutrondetector.

DETAILED DESCRIPTION OF SPECIFIC EMBODIMENTS

One or more specific embodiments of the presently disclosed subjectmatter are described below. In an effort to provide a concisedescription of these embodiments, not all features of an actualimplementation are described in the specification. It should beappreciated that in the development of any such actual implementation,as in any engineering or design project, numerousimplementation-specific decisions must be made to achieve thedevelopers' specific goals, such as compliance with system-related andbusiness-related constraints, which may vary from one implementation toanother. Moreover, it should be appreciated that such a developmenteffort might be complex and time consuming, but would nevertheless be aroutine undertaking of design, fabrication, and manufacture for those ofordinary skill having the benefit of this disclosure.

Embodiments of the presently disclosed subject matter generally relateto systems and methods for well logging using gamma-rays produced byinelastic scattering events (“inelastic gamma-rays”). Counts or spectraof inelastic gamma-rays may indicate a variety of properties of asurrounding subterranean formation. For example, in combination with ahydrogen index signal, a count of inelastic gamma-rays may enabledetermination of porosity and gas saturation. FIG. 1 illustrates asystem 10 for determining formation properties using inelasticgamma-rays that includes a downhole tool 12 and a data processing system14. By way of example, the downhole tool 12 may be a slickline orwireline tool for logging an existing well, or may be installed in aborehole assembly for logging while drilling (LWD). The data processingsystem 14 may be incorporated into the downhole tool 12 or may be at aremote location.

The downhole tool 12 may be constructed so as to improve detection ofgamma-rays produced via inelastic scattering events, while reducingdetection of gamma-rays produced via thermal and epithermal neutroncapture events. Indeed, the downhole tool 12 may provide for a gamma-rayresponse substantially free of epithermal and thermal neutron capturebackground, such that the gamma-ray response is substantiallyindependent of neutron porosity. With reference to FIG. 2, the totalgamma-ray response 1 of an example downhole tool is decomposed by MonteCarlo modeling into the inelastic scattering 2, epithermal capture 3,and thermal capture 4 components. For the particular components in thistool, the sum of epithermal capture counts during the burst isapproximately ⅓ of the sum of inelastic counts and is, therefore, quitelarge. While it might appear that one could estimate the epithermalcapture amplitude and decay rate by sampling the count rate immediatelyafter the burst and with that, correct for the epithermal componentduring the burst, this is difficult in practice because the epithermalcapture count rate after the burst is much smaller than the thermalcapture count rate. In effect, the thermal capture count rate during theburst-off period overwhelms the epithermal capture count rate, such thatthe epithermal count rate and decay rate cannot be determinedindependently. The fraction of the epithermal capture component relativeto the inelastic scattering component varies depending on the materialsused in the downhole tool The particular materials employed in each ofthe components may be described in greater detail below, with particularreference to FIG. 3. With continued reference to FIG. 1, the downholetool 12 may be surrounded by a housing 16 made of elements with lowresonance integrals such as cobalt-free steel. A neutron source 18 maybe any suitable neutron source capable of emitting neutrons into asurrounding formation to produce inelastic gamma-rays. By way ofexample, the neutron source 18 may be a pulsed electronic neutronsource, such as a Minitron™ by Schlumberger Technology Corporation.Additionally or alternatively, in certain embodiments, the neutronsource 18 may be a radioisotope source capable of emitting fastneutrons. A neutron shield 20 containing elements with high (n,2n) crosssections, such as lead, bismuth or tungsten or combinations of these,may separate the neutron source 18 from various detectors in thedownhole tool 12.

Certain embodiments of the downhole tool 12 may include a neutrondetector 22, which may be surrounded by a cadmium-containing thermalneutron shield 24. The neutron detector 22, which may be a ³He neutrondetector, may primarily detect epithermal neutrons rather than thermalneutrons, since the thermal neutron shield 24 may serve to prevent thepassage of thermal neutrons to the detector 22. The downhole tool 12 mayoptionally include a neutron monitor 18.1, which is located near theMinitron and detects primarily unscattered neutrons directly from theneutron generator. The optional neutron monitor, which may be a plasticscintillator coupled to a photomultiplier, provides a count rate signalproportional to the neutron output rate from the generator. The downholetool 12 may include one or more gamma-ray detectors, illustrated in FIG.1 as a “near” gamma-ray detector 26 and a “far” gamma-ray detector 28,each of which may be surrounded by a housing 30. As should beappreciated, the near gamma-ray detector 26 and the far gamma-raydetector 28 are so named due to their relative proximity to the neutronsource 18. In some embodiments, a scintillation crystal 32 of the neargamma-ray detector 26 may be located approximately 5″ to 22″ from theneutron source 18, while the scintillation crystal 32 of the fargamma-ray detector 28 may be located approximately 15″ to 36″ from theneutron source 18.

As noted above, the near gamma-ray detector 26 and the far gamma-raydetector 28 may be contained in respective housings 30. For reasonsdescribed below, each of the housings 30 may contain materialsincorporating ⁶Li, such as lithium carbonate (Li₂CO₃), which maysubstantially shield the gamma-ray detectors 26 and 28 from thermalneutrons without producing thermal neutron capture gamma-rays. Thescintillation crystals 32 of the gamma-ray detectors 26 and 28 mayenable detection counts or spectra of gamma-rays by producing light whengamma-rays are scattered or absorbed by the scintillation crystals 32.For reasons described below, housings 34 of aluminum alloy or fiberglassmay surround the scintillation crystals 32 to reduce production ofepithermal neutron capture gamma-rays. Depending on the application, thescintillation crystals 32 may be chosen to include any of a variety ofmaterials, as described below. Photodetectors 36 may detect lightemitted by the scintillation crystals 32 when a gamma-ray is absorbed,once the light has passed through an optical window 38, to obtain agamma-ray count.

The signals from the neutron detector 22, the near gamma-ray detector26, and/or the far gamma-ray detector 28 may be transmitted to the dataprocessing system 14 as data 40. The data processing system 14 mayinclude a general-purpose computer, such as a personal computer,configured to run a variety of software, including software implementingall or part of the present technique. Alternatively, the data processingsystem 14 may include, among other things, a mainframe computer, adistributed computing system, or an application-specific computer orworkstation configured to implement all or part of the present techniquebased on specialized software and/or hardware provided as part of thesystem. Further, the data processing system 14 may include either asingle processor or a plurality of processors to facilitateimplementation of the presently disclosed functionality.

In general, the data processing system 14 may include data processingcircuitry 44, which may be a microcontroller or microprocessor, such asa central processing unit (CPU), which may execute various routines andprocessing functions. For example, the data processing circuitry 44 mayexecute various operating system instructions as well as softwareroutines configured to effect certain processes and stored in orprovided by a manufacture including a computer readable-medium, such asa memory device (e.g., a random access memory (RAM) of a personalcomputer) or one or more mass storage devices (e.g., an internal orexternal hard drive, a solid-state storage device, CD-ROM, DVD, or otherstorage device). In addition, the data processing circuitry 44 mayprocess data provided as inputs for various routines or softwareprograms, including the data 40.

Such data associated with the present techniques may be stored in, orprovided by, the memory or mass storage device of the data processingsystem 14. Alternatively, such data may be provided to the dataprocessing circuitry 44 of the data processing system 14 via one or moreinput devices. In one embodiment, data acquisition circuitry 42 mayrepresent one such input device; however, the input devices may alsoinclude manual input devices, such as a keyboard, a mouse, or the like.In addition, the input devices may include a network device, such as awired or wireless Ethernet card, a wireless network adapter, or any ofvarious ports or devices configured to facilitate communication withother devices via any suitable communications network, such as a localarea network or the Internet. Through such a network device, the dataprocessing system 14 may exchange data and communicate with othernetworked electronic systems, whether proximate to or remote from thesystem. The network may include various components that facilitatecommunication, including switches, routers, servers or other computers,network adapters, communications cables, and so forth.

The downhole tool 12 may transmit the data 40 to the data acquisitioncircuitry 42 of the data processing system 14 via, for example, atelemetry system communication downlink or a communication cable. Afterreceiving the data 40, the data acquisition circuitry 42 may transmitthe data 40 to data processing circuitry 44. In accordance with one ormore stored routines, the data processing circuitry 44 may process thedata 40 to ascertain one or more properties of a subterranean formationsurrounding the downhole tool 12. Such processing may involve, forexample, one or more techniques for reducing an epithermal neutroncapture background from a gamma-ray count, as described below. The dataprocessing circuitry 44 may thereafter output a report 46 indicating theone or more ascertained properties of the formation, such as porosityand gas saturation, as discussed below. The report 46 may be stored inmemory or may be provided to an operator via one or more output devices,such as an electronic display and/or a printer.

FIG. 3 illustrates a well-logging operation 48 involving the downholetool 12 in a surrounding subterranean formation 50. In the operation 48depicted in FIG. 3, the downhole tool 12 has been lowered into anexisting well surrounded by casing 52. The well-logging operation 48 maybegin when the neutron source 18 outputs a burst of neutrons 54 ofapproximately 14.1 MeV or greater into the surrounding formation 50.Initially, the burst of neutrons 54 may collide with nuclei of theformation 50 in inelastic scattering events 56, which causes inelasticgamma-rays 58 to be emitted and the neutrons of the burst of neutrons 54to lose energy. Following or concurrently with the inelastic scatteringevents 56, the neutrons of the burst of neutrons 54 may interact withthe formation 50 in elastic scattering events 60, which causes theneutrons to drop in energy to eV and sub-eV levels as epithermalneutrons 62 and/or thermal neutrons 64.

The inelastic gamma-rays 58 produced by the inelastic scattering events56 may subsequently Compton-scatter due to electrons in the formation50. Some of the inelastic gamma-rays 58 may ultimately be detected inthe near gamma-ray detector 26 or the far gamma-ray detector 28. Theresulting count of the gamma-rays 58 may be sensitive to formationdensity for at least two reasons. First, the probability that one of theneutrons of the burst of neutrons 54 will interact with a nucleus of theformation 50 may be proportional to the number density of nuclei in theformation 50. Thus, if the formation 50 includes a greater numberdensity of nuclei for the same chemical composition of the formation,the mean path between neutron collisions with nuclei becomes shorter.The neutrons lose energy through collisions faster and travel a shorteraverage distance from the neutron source before dropping below 1 MeV.Once the neutrons drop below approximately 1 MeV, the neutrons no longerhave sufficient energy to create inelastic gamma-rays 58. Thus, agreater number density in the formation will decrease the number ofinelastic gamma-rays produced far from the neutron source and fewerinelastic gamma-rays will be produced near the far detector 28. Sincegamma-rays produced near the far detector 28 have a high likelihood ofbeing detected, the far detector count rate is very sensitive to thisneutron transport effect. Because the near detector 26 is closer to theneutron source, the decrease in the number of inelastic gamma-raysproduced near the detector is much smaller than for the far detector 28.Second, the inelastic gamma-rays 58 created by the inelastic scatteringevents 56 may Compton-scatter more frequently if the formation 50 ismore dense. As a result, the gamma-rays 58 may be less likely to reachthe near gamma-ray detector 26 or the far gamma-ray detector 28 beforebeing photoelectrically absorbed by the formation 50.

Either of the above-described effects may lead to a lower inelasticgamma-ray count as the density of the formation 50 increases. It shouldbe noted, however, that the former effect is sensitive to the nucleusnumber density of the formation 50, while the latter effect is sensitiveto the electron number density of the formation 50 or, approximately,the bulk density of the formation 50. The nucleus number density may notcorrelate well with bulk density because heavy nuclei and light nucleioften have similar neutron interaction probabilities; however, thecombination of nuclei number density and bulk density obtained from aneutron-gamma measurement may be used in much the same way as aconventional density measurement based on Compton scattering ofgamma-rays.

A complication may arise, however, due to the interaction of theneutrons of the neutron burst 54 with hydrogen in the formation 50.Because elastic scattering 60 occurs primarily due to interactions withhydrogen in the formation 50, the number of low-energy epithermalneutrons 62 and thermal neutrons 64 that reach the vicinity of thegamma-ray detectors 26 and 28 may accordingly be strongly influenced bythe hydrogen index of the formation 50. If one of the epithermalneutrons 62 or thermal neutrons 64 were captured by a nucleus at or inthe vicinity of the near gamma-ray detector 26 or the far gamma-raydetector 28, such as in the casing 52, a neutron capture event 66 mayoccur. Such a neutron capture event 66 may produce a neutron capturegamma-ray 68 that may be detected by the gamma-ray detector 26 or 28.Without correction, neutron capture gamma-rays 68 may completelyoverwhelm the inelastic gamma-rays 58, resulting in gamma-ray countswith the character of neutron porosity measurements, rather than densitymeasurements. Such measurements would be of less use to combine with aneutron porosity (hydrogen index) measurement to identify gas becausethe two measurements would be highly correlated through hydrogen index.

For this reason, the downhole tool 12 may contain materials carefullychosen to reduce or eliminate neutron capture events 66 occurring in thedownhole tool 12. Moreover, because some neutron capture events 66 mayyet still occur in the downhole tool 12, and because some neutroncapture events 66 may take place externally to the downhole tool 12,various techniques may be employed to reduce or eliminate any remainingneutron capture background from gamma-ray signals of the gamma-raydetectors 26 and/or 28. In particular, the materials chosen and thetechniques employed may reduce or eliminate not only thermal neutroncapture background, but also epithermal neutron capture background, fromgamma-ray signals detected by the gamma-ray detectors 26 and/or 28.Additionally, certain materials in the downhole tool 12, such as thescintillation crystals 32, may be chosen based on additional criteriaparticular to a chosen well-logging application.

Considering first the material chosen for the scintillation crystals 32in the gamma-ray detectors 26 and/or 28, the capture probability ofepithermal neutrons in matter may be described approximately by a valueknown as “resonance integral,” which is the integral cross section forinteraction in the neutron resonance energy (“epithermal”) region. Theresonance integral for detector materials of interest are shown below inTable 1 (from “Nuclides and Isotopes Chart of the Nuclides,” KnollsAtomic Power Laboratory, 16^(th) Edition). The thermal capture crosssection is indicated by the parameter σ_(a). The materials described inTable 1 may generally form a component of a scintillation crystal.

TABLE 1 σ_(a) Resonance Integral Element (barns) (barns) RelevantMaterial Gd 49000 400 GSO La 9 12 Lanthanum halide detectors Si 0.1680.1 GSO, Glass scintillators O 0.0 0.0 GSO, Glass scintillators Na 0.5332 NaI detector Cl 33.5 12 LaCl₃ detector Br 6.8 89 LaBr₃ detector I 6.2150 NaI Lu 84 900 Lu-based detectors Y 1.28 1 Y-based detectors (e.g.,YAlO₃) Al 0.23 0.17 YAP detectors Bi 0.034 0.19 BGO detectors Ge 2.2 6BGO detectors Ce 0.63 0.7 Activator in LaCl₃, LaBr₃, GSO, YAP, LuAP

From Table 1, it may be seen that certain materials all have componentswith a “resonance integral” of 89 or larger, such as GSO (gadoliniumoxyorthosilicate, or Gd₂SiO₅); lutetium-based scintillators such as LuAP(LuAlO₃); LaBr₃; and NaI. On the other hand, certain other materialshave no component with a resonance integral greater than 12, making themlargely insensitive to epithermal neutron capture. Such materialsinclude BGO (bismuth germanate, or Bi₄Ge₃O₁₂); LaCl₃; and YAP (YAIO₃).

The significance of the resonance integral and thermal capture crosssection will be appreciated since neutron capture within thescintillation crystal generally results in the emission of a gamma-ray,which then has a high probability of being absorbed and detected by thescintillator. Depending on the particular nuclei in the scintillator andin the surrounding tool and formation, different energies of gamma-rayswill be absorbed within the scintillator. In some applications whereonly high energy gamma-rays are desired, it is common to discriminateagainst lower energy gamma-rays with a pulse-height discriminator.However, in the application of gas detection and measurement, it isdesirable to detect all inelastic gamma-rays entering the scintillator,regardless of energy, to optimize precision. It is, therefore, notdesirable to reduce neutron capture gamma-rays by means of pulse heightdiscrimination as this would also eliminate some inelastic gamma-rays.As a result, minimizing the resonance integral is meaningful.

In choosing the material for the scintillation crystals 32, otherfactors may be considered. Specifically, certain applications, such asCarbon/Oxygen (C/O) logging, may involve detection of inelasticgamma-rays 58 for use in spectroscopy analyses (“inelasticspectroscopy”). For such applications, the scintillation crystals 32 maybe chosen to obtain an optimal precision (relative logging speed) forinelastic spectroscopy that may also optimize the response for the moststatistically significant number of elemental yields. In an extensivestudy of scintillation crystals of various materials, the results ofwhich are provided in Table 2 below, several factors pertinent toprecision (relative logging speed) were considered according to thefollowing relationship:

(Relative Logging Speed)˜(Relative Spectral Quality)*(RelativeEfficiency)*(Relative Maximum Counting Rate)*(Relative Improvement Dueto Neutron Background Removal)  (1).

As considered in Equation (1) above, Relative Spectral Quality (highervalue is better) represents a measure of the ability to separate, by aleast-squares process, the spectral signatures of different elements ina statistical manner. Relative Spectral Quality includes manyscintillation crystal properties, such as light output, atomic number,temperature response, peak-to-Compton ratio, size, and resolution.Relative Efficiency (higher value is better) represents a measure of thefraction of high-energy gamma-rays that interact while passing throughthe scintillation crystal and that are therefore detected. This valuerelates to scintillation crystal size, density, and atomic number.Relative Maximum Counting Rate (higher value is better) represents ameasure of how fast the scintillation crystal is able to detect andprocess individual gamma-rays that are absorbed in the scintillator, andis based on the light production and decay properties of thescintillator. Relative Improvement due to Neutron Background Removal(higher value is better) is a measure of the statistical improvement inthe measurement due to a change in the background neutron signal in thescintillation crystal. This number will increase with decreasing neutronbackground.

TABLE 2 Relative Relative Relative Relative Relative ImprovementRelative Logging Speed Logging Speed Spectral High-Energy due to NeutronMax (no neutron (neutron Detector Quality Efficiency Bkg RemovalCountrate limit) limited) LaCl 1.63 0.32 3.01 11.5 12.5 1.7 LaBr 1.620.39 1.73 6.6 5.0 1.2 NaI 1.65 0.32 1.73 1.0 1.0 1.0 GSO 1.84 0.66 1.003.9 4.9 1.3 BGO 2.60 0.61 2.33 0.77 2.3 4.0 LuAP 1.39 0.65 1.03 13.612.1 1.0

Table 2 illustrates how these properties may vary for scintillationcrystals of the same size of various types. The last two columns ofTable 2 describe the relative logging speed (higher value is better) forthe detector when all other factors are the same. Specifically, relativelogging speed (no neutron limit) representing when measurements takenwith a scintillation crystal are not neutron limited, meaning thatenough neutrons are produced so as to push the scintillation crystal toits limit. Relative logging speed (neutron limited) provides a valuerepresenting when measurements taken with a scintillator are neutronlimited, meaning that fewer neutrons are produced than a maximumcapability of the scintillator. The numbers provided in Table 2 arebased on a least-squares processing method, but it should be appreciatedthat the data may also be processed using a standard “windows”processing.

It should be appreciated that, in optimizing a tool containing one ofthe scintillation crystals listed above in Table 2, one would alsoadjust the position of the scintillation crystal relative to the neutronsource, which may optimize the countrate of the scintillation crystalversus the degradation in formation response. Therefore, an optimum tooldesign may have an effective relative logging speed somewhere betweenthe values in the two columns. It should also be noted that the spectralquality factor, neutron background term, and efficiency may change asthe size of the scintillation crystal changes, which may also affect thevalues in the last two columns. As apparent in Table 2, spectral qualitymay not significantly impact the values of relative logging speed, butthe three remaining factors may significantly impact the values ofrelative logging speed.

Based on the above considerations, the scintillation crystal 32 in thegamma-ray detectors 26 and/or 28 of the downhole tool 12 may includematerials with a resonance integral less than 12, such as BGO (bismuthgermanate, or Bi₄Ge₃O₁₂); LaCl₃; and YAP (YAlO₃). In particular, thescintillation crystals 32 may be formed of a lanthanum-halidescintillation crystal, such as LaCl, for inelastic gamma-rayspectroscopy applications, or of YAP (YAlO₃), for applications relianton inelastic gamma-ray counts, such as porosity and gas saturationmeasurements. With such materials, epithermal neutrons 62 that passthrough components of the downhole tool 12 to reach one of thescintillation crystals 32 may be less likely to cause neutron captureevents 66 within the scintillation crystals 32. Accordingly, epithermalneutron capture gamma-ray background may be reduced.

By choosing the scintillation crystal 32 such that the epithermalneutron capture background is low enough such that the epithermalcomponent does not have to be removed by other means, an inelasticgamma-ray spectroscopy measurement may be improved due to variousfactors. First, if the epithermal capture background is not removed, thestatistical variance (error) may not be inflated due to the process ofsubtracting the epithermal neutron capture background. Second, removalof unwanted counts may enable the designer of the downhole tool 12 toincrease neutron output from the neutron source 18 or to change thespacing of the near and far gamma-ray detectors 26 and/or 28 such thatthe maximum count rate capability of each can be utilized for usefulcounts. Third, some potential complicated neutron physics may be removedfrom the petrophysical interpretation. The resulting better performancemay be indicated by improved statistical precision and, hence, fasterlogging speed, by the ability to extract more statistically validinelastic elements, and by the reduction of unwanted physics from theanalysis. This improves not only the statistical precision but also theaccuracy of the measurement. To reduce the effect of neutron capture inthe gamma-ray detectors 26 and/or 28, the scintillation crystals 32 maybe surrounded by thermal neutron shields 34. Specifically, to reducethermal neutron capture gamma-ray production, the thermal neutronshields 34 may employ a material enriched with the isotope ⁶Li, whichhas a relatively high thermal neutron capture cross section, but whichalso produces primarily charged particles in lieu of gamma-rays as aconsequence of the capture reaction ⁶Li(n,α)³H. Candidate materialsinclude lithium carbonate (Li₂CO₃) embedded in epoxy, or metalliclithium or Li₃N hermetically sealed in a metal can. As illustrated inFIG. 4, when the thermal neutron shields 34 include such a material,thermal neutrons 64 that reach the thermal neutron shields 34 may beabsorbed without the production of neutron capture gamma-rays 68. OtherLi compounds could be considered including Li₃PO₄, TetraLithium EDTA orless desirable LiF. An appropriate Li alloy may also be considered if itis formulated to reduce the natural chemical reactivity of the containedLi.

Turning briefly to FIG. 4, a plot 79 illustrates the neutron shieldingcapabilities of lithium carbonate (Li₂CO₃) enriched with the isotope⁶Li. In the plot 79, an ordinate 80 represents a fraction of epithermalneutrons 62 and thermal neutrons 64 passing through a quantity oflithium carbonate (Li₂CO₃), displayed logarithmically from 10⁻² to 10⁰.An abscissa 81 represents the thickness of the lithium carbonate(Li₂CO₃) in units of millimeters (mm). Data points 82 and 84 representquantities of thermal neutrons 64 and epithermal neutrons 62,respectively, that pass through the lithium carbonate (Li₂CO₃). As shownin the plot 79, when the thickness of the lithium carbonate (Li₂CO₃)increases, the fraction of thermal neutrons 64 passing through thelithium carbonate (Li₂CO₃) drops dramatically. Accordingly, when thethermal neutron shields 34 include lithium carbonate (Li₂CO₃) enrichedwith the isotope ⁶Li, the quantity of thermal neutrons 64 that reach thescintillation crystals 32 may be substantially reduced.

An additional embodiment would be to add a third detector to thedownhole tool of FIG. 1 for spectroscopy purposes, such that it would beplaced between detectors 26 and 28. The scintillator of this detectorcould be chosen from materials of higher resonance integral as long asgamma-rays from capture of epithermal neutrons within the scintillatorwere not detected in significant amounts in detectors 26 and 28. Toprevent such gamma-rays from being detected in detectors 26 and 28,shields made of lead, bismuth, tungsten or a combination of the threemay be placed between the third detector and each of detectors 26 and28.

Because epithermal neutrons 62 may also reach other components of thedownhole tool 12, such as the housing 16 or the neutron shield 20, othercomponents of the downhole tool 12 may be constructed of materialshaving a low resonance integral. The resonance integral for certainrelevant materials are listed below in Table 3.

TABLE 3 σ_(a) Resonance Integral Element (barns) (barns) RelevantMaterial Fe 2.56 1.4 Steel Ni 4.5 2.2 Steel Cr 3.1 1.6 Steel Co 37.2 74Steel Nb 1.15 8.5 Steel Mo 2.5 25 Steel Mn 13.3 14 Steel V 5 2.8 SteelTi 6.1 2.9 titanium Al .23 .17 Aluminum alloy Mg .066 .038 Aluminumalloy C 0 0 Fiberglass O 0 0 Fiberglass H .332 .149 Fiberglass Si .168.1 Fiberglass W 18.2 350 Shield Pb .171 .14 Shield Bi .034 .19 Shield

As shown in Table 3, the resonance integral of steels may be relativelylow. However, the production of neutron capture gamma-rays 68 fromneutron capture events 66 in steel can be high in large quantities ofsteel. Thus, where possible, the housing 16 and other various componentsof the downhole tool 12 should be made of aluminum alloy or fiberglass.To the extent that steel is used in components of the downhole tool 12,the steel should contain very little or no cobalt (Co), manganese (Mn),molybdenum (Mo) or niobium (Nb) to reduce the occurrence of epithermalneutron capture in the downhole tool 12. Thus, with reference to FIG. 3,the housing 16 of the downhole tool 12 may be constructed of steelwithout the previously mentioned elements, and the housings 30 of thegamma-ray detectors 26 and/or 28 may be constructed primarily ofaluminum alloy or fiberglass. The thermal neutron shield 24 that maysurround the neutron detector 22 may contain cadmium (Cd), and theneutron shield 20 between the neutron source 18 and the various othercomponents of the downhole tool 12 may contain tungsten (W), lead (Pb),)or bismuth (Bi). Tungsten has desirable properties as both a gamma-rayshield and a neutron shield but it has a high resonance integral. Toprevent gamma-rays from epithermal neutron capture in a tungsten shieldfrom entering the detectors 26 and 28, a shield of lead (Pb) or bismuth(Bi) may be placed between the tungsten and detectors to shield thedetectors from epithermal neutron capture gamma-rays from tungsten.

For some choices of scintillators 32, it may be desirable to includegamma-ray shielding between detectors 26 and 28. For example, if a LaCl₃scintillator is chosen for detector 26 and a YAP detector is chosen fordetector 28, some gamma-rays from neutron capture in the Cl of the LaCl₃detector could be detected in the YAP scintillator. To prevent this,gamma-ray shielding composed of lead (Pb), bismuth (Bi), or tungsten (W)can be placed between the detectors. Additionally, if the neutrondetector 24 has a thermal neutron filter that emits gamma-rays (such ascadmium (Cd) or Gadolinium (Gd)), a gamma-ray shield composed of Pb, Bi,or W can be placed between detectors 24 and 26 to prevent capturegamma-rays from the neutron filter from entering detector 26. To preventgamma-rays from epithermal neutron capture in a tungsten shield fromentering the detectors 26 and 28, a shield of lead (Pb) or bismuth (Bi)may be placed between the tungsten and detectors to shield the detectorsfrom epithermal neutron capture gamma-rays from tungsten.

With the components of the downhole tool 12 optimized to reduce thenumber epithermal neutron capture events 66 in the downhole tool 12, thedownhole tool 12 may be employed for a variety of applications. Forexample, FIG. 5 provides a flowchart 86 describing an embodiment of amethod for logging porosity and gas saturation in a subterraneanformation, such as the subterranean formation 50 of FIG. 3. In a firststep 88, a burst of neutrons 54 of approximately 14.1 MeV may be emittedinto the formation 50 via the neutron source 18. The emitted burst ofneutrons 54 may interact with nuclei of the formation 50, resulting ininelastic scattering events 56 that produce inelastic gamma-rays 58. Theburst of neutrons 54 may continue to lose energy through elasticscattering events 60, losing energy to become epithermal neutrons 62 andthermal neutrons 64.

In step 90, epithermal neutrons 62 that scatter back toward the downholetool 12 may be detected using the neutron detector 22 surrounded by thethermal neutron shield 24 to obtain an epithermal neutron count. Asnoted above, the thermal neutron shield 24 may exclude substantially allthermal neutrons 64 from entering the neutron detector 22, resulting ina neutron signal that includes a count of substantially only epithermalneutrons 62. Alternatively, the epithermal neutron count may be obtainedby detecting both epithermal neutrons 62 and thermal neutrons 64 usingthe neutron detector 22 without the thermal neutron shield 24 to obtaina total neutron count, before correcting the total neutron count toinclude only epithermal neutrons, as described below with reference toFIG. 8.

In step 92, the near gamma-ray detector 26 and/or the far gamma-raydetector 28 may detect inelastic gamma-rays 58 from the formation 50, aswell as neutron capture gamma-rays 68 that may still occur in or aroundthe downhole tool 12, to obtain a raw gamma-ray count. Though the rawgamma-ray count may consist substantially of inelastic gamma-rays 58, aneutron capture background due to neutron capture gamma-rays 68 maycontaminate the signal. FIG. 6 shows two burst periods of the totalgamma-ray multichannel scalar count rate 1, inelastic gamma-ray countrate 2, the epithermal capture gamma-ray count rate 3, the thermalcapture gamma-ray count rate 4, and the sum of epithermal and capturegamma-ray count rate (total capture gamma-ray count rate) 5 for adownhole pulsed-neutron tool when low-resonance-integral materials arechosen according to this invention. The epithermal capture gamma-raybackground is greatly reduced compared to that in FIG. 1.1 of the priorart downhole tool. However, a small epithermal capture gamma-raybackground and a significant thermal capture gamma-ray backgroundremains. Thus, in step 94, the raw gamma-ray count may be corrected toreduce or eliminate the thermal and epithermal neutron capturebackground, such that the resulting gamma-ray count includessubstantially only inelastic gamma-rays 58. The correction of step 94may involve any of a variety of techniques, many of which are describedbelow.

The total gamma-ray count rate illustrated in FIG. 6 has a sharp risetime at the beginning of a burst (A) and a sharp fall time at the end ofthe burst (B). Times A and B are straightforward to determine by variousmeans; for example, from the sharp breaks in the count rate with time orby fitting regions B and C as defined in FIG. 7. Between points A and B,the thermal capture gamma-ray count rate 4 first decreases from itsvalue at A and then increases as time approaches B. It continues toincrease after the end of the burst for a short time as the neutronsemitted during the burst slow down to thermal energy. The epithermalcount rate 3 has a complementary behavior. It first increases from itsvalue at A (zero) and then becomes constant as time approaches B.Finally it decreases after the end of the burst for a short time as thelast of the neutrons emitted during the burst decrease in energy tothermal. Because of the complementary nature of the epithermal 3 andthermal 4 capture gamma-rays during the burst, the sum of the twocomponents 5 is almost linear between points A and B. Referring to FIG.10, a simple linear connection 2 of the count rate at A to B is a closerepresentation of the combined epithermal and thermal capture gamma-raycount rate 3. The sum of counts under the linear connection 2 is 14%lower than the sum of the total capture counts. Therefore, an accuraterepresentation of the sum of the total capture gamma-ray counts is 1.14times the sum of the linear connection counts. In practice the actualrepresentation of the sum of the total capture gamma-ray counts is ftimes the sum of the linear connection counts, where f is a number maybe different for each detector type and spacing from the source and isdetermined by calibration. Once the sum of the total capture gamma-raycounts is determined, it may be subtracted from the total gamma-raycounts during the burst to yield the sum of inelastic gamma-rays duringthe burst. An alternative to a linear background subtraction is to use apolynomial background subtraction where the additional terms in thepolynomial are determine by requiring that the inelastic count rateafter background subtraction is constant or nearly constant during theburst, or by Monte Carlo modeling.

It should be appreciated that the previous method is best applied insituations where the epithermal capture gamma-ray background is smallcompared to the inelastic gamma-ray component as a consequence of thelow-resonance-integral nature of this invention. However, in thefollowing several other capture gamma-ray subtraction methods aredescribed as alternative embodiments that may be appropriate even if theepithermal capture background is not small compared to the inelasticgamma-ray component.

Having obtained a signal which may be a count of epithermal neutrons ormay be a ratio of capture gamma-rays between the near 26 and far 28detectors, which corresponds to neutron porosity and to hydrogen index,and a count of inelastic gamma-rays 58, which may correspond to aformation density measurement, the data processing circuitry 14 maydetermine porosity and gas saturation in step 96. One manner ofdetermining porosity and gas saturation using the count of inelasticgamma-rays 58 is described below with particular reference to FIG. 13

FIG. 7 is a flowchart 98 representing an embodiment of a method forreducing or eliminating an epithermal neutron capture background fromthe raw gamma-ray count. The flowchart 98 may generally represent onemanner of performing step 94 of the flowchart 86 of FIG. 7. In a firststep 100, the data processing circuitry 44 of the data processing system14 may receive a count of epithermal neutrons 62, which may be obtainedfrom the neutron detector 22 surrounded by the thermal neutron shield24. As described above, because the thermal neutron shield 24 largelyprevents the detection of thermal neutrons 64, while permitting thedetection of epithermal neutrons 62, the neutron count may includesubstantially only epithermal neutrons 62 without requiring furtherprocessing.

In step 100, the data processing circuitry 44 of the data processingsystem 14 may obtain a count of inelastic gamma-rays 58 by reducing theraw gamma-ray count by an amount proportional to the epithermal neutroncount of step 98. Thus, the corrected gamma-ray count may includesubstantially only inelastic gamma-rays 58, and may therefore beindependent of the neutron porosity measurement.

FIG. 8 is a flowchart 104 representing an embodiment of an alternativemethod for reducing or eliminating an epithermal neutron capturebackground from the raw gamma-ray count. The flowchart 104 may alsogenerally represent one manner of performing step 94 of the flowchart 86of FIG. 7. In a first step 106, the data processing circuitry 44 of thedata processing system 14 may receive a count of both epithermalneutrons 62 and thermal neutrons 64, which may be obtained from aneutron detector 22 that lacks a thermal neutron shield 24. The neutroncount ascertained in step 106 may represent all neutrons detected at atime when the neutron source 18 is emitting the neutron burst 54. Duringthis time, the short-lived epithermal neutrons 62 may be presentalongside the longer-lived thermal neutrons 64.

In step 108, the data processing circuitry 44 may receive a count ofthermal neutrons, which may represent all neutrons detected by theunshielded neutron detector 22 between emissions of the neutron burst54. Because the thermal neutrons 64 are longer-lived than the epithermalneutrons 62, substantially all of the neutrons detected after theneutron burst 54 may be thermal neutrons. Thus, the neutron countreceived in step 108 may accordingly include substantially only thermalneutrons 64.

In step 110, the data processing circuitry 44 may determine anapproximate epithermal neutron count by reducing the neutron count ofstep 106 by the normalized neutron count of step 108. In step 112, thedata processing circuitry 44 may determine a count of inelasticgamma-rays 58 by reducing the raw gamma-ray count by an amountproportional to the approximate epithermal neutron count of step 110.The corrected gamma-ray count may therefore include substantially onlyinelastic gamma-rays 58.

Alternative methods of performing the correction of step 94 of FIG. 5may not involve a count of epithermal neutrons 62, as described abovewith reference to FIGS. 5 and 6, but may rather rely upon the temporalnature of the gamma-ray response during and following the burst ofneutrons 54. Turning to FIG. 9, a plot 114 illustrates a simulatedgamma-ray count from a single gamma-ray detector, such as the neargamma-ray detector 26, during and after the burst of neutrons 54. Anordinate 116 represents the total count weight of the detectedgamma-rays and an abscissa 118 represents time in units of microseconds(μs). A curve 120 represents the total count of gamma-rays at variouspoints in time.

A first time segment between 0 and 20 microseconds, labeled A,represents a time when the neutron burst 54 is being emitted from theneutron source 18. A second time segment from 20 microseconds to 40microseconds, labeled B, represents a time following the neutron burst54. In the time segment B, nearly all inelastic scattering events 56have ended, and the detected gamma-rays are due substantially entirelyto capture events 66 involving epithermal neutrons 62 and thermalneutrons 64. A third time segment between 40 microseconds and 100microseconds, labeled C, represents a time when substantially allgamma-rays detected are due to capture events 66 involving thermalneutrons 64. Based on data represented by the curve 120, the dataprocessing circuitry 44 of the data processing system 14 may reduce oreliminate epithermal neutron capture background from a raw gamma-raycount obtained by the near gamma-ray detector 26 and/or the fargamma-ray detector 28. Such techniques may be described with referenceto FIGS. 11 and 12.

Turning to FIG. 11, a flowchart 122 represents another embodiment of amethod for correcting for, reducing, or eliminating an epithermalneutron capture background from a raw gamma-ray count. The flowchart 122may generally represent one manner of performing step 94 of theflowchart 86 of FIG. 5. In a first step 124, the data processingcircuitry 44 of the data processing system 14 may receive a count ofgamma-rays detected for time segments A, B, and C. Such a gamma-raycount may represent both inelastic gamma-rays 58 as well as epithermaland thermal neutron capture gamma-rays 68.

In step 126, the data processing circuitry 44 may fit an exponentialfunction to the curve 120 in the intervals B and C. The exponentialfunction may include a term for thermal neutron capture, fitted in theinterval C, and an additional term for epithermal neutron capture,fitted in the interval B. In this way, terms proportional to theepithermal count rate and thermal count rate may be determinedseparately. In step 128, the data processing circuitry 44 may apply theexponential function to the interval A. The exponential function may beused for subtracting appropriate fractions of the epithermal count rateand the thermal count rate from the gross count rate in the interval A,leaving a corrected count of substantially only inelastic gamma-rays 58.

In lieu of, or in addition to, fitting an exponential function to thecurve 120 to ascertain thermal and epithermal neutron capturebackground, the data processing circuitry 44 may employ inverse Laplacetransforms to ascertain the same. Such techniques are described ingreater detail with regard to compensation for thermal neutron capturebackground in U.S. application Ser. No. 11/768,034, “Well Logging MethodFor Determining Formation Characteristics Using Pulsed Neutron CaptureMeasurements,” (Atty. Docket No. 60.1672) which is incorporated byreference herein in its entirety. The techniques may be extended tocompensate for epithermal neutron capture background by considering, asdescribed by the flowchart 122 of FIG. 11, that the interval Bimmediately following the neutron burst 54 includes both epithermal andthermal neutron capture background, while the subsequent interval Cincludes substantially only thermal neutron capture background.

FIG. 12 is a flowchart 130 representing another embodiment of a methodfor reducing or eliminating an epithermal neutron capture backgroundfrom a raw gamma-ray count. The flowchart 130 may generally representone manner of performing step 94 of the flowchart 86 of FIG. 5 In theflowchart 130, epithermal neutron capture background is corrected basedon the increase in gamma-ray count that occurs during the interval A.

In a first step 132, the data processing circuitry 44 may receive acount of the gamma-rays detected during the interval A, which representsthe time interval during which the neutron burst 54 is emitted. In step134, the data processing circuitry 44 may determine the quantity ofepithermal and thermal neutron capture gamma-rays 68 by analyzing theincrease in raw gamma-ray count over the interval A. It should beappreciated that if the neutron output during the burst 54 is notconstant as a function of time, a correction can be made using aninstantaneous count rate from a fast neutron monitor that may be placednear the neutron source 18. Turning briefly to FIG. 9, the curve 120indicates a sharp rise at the beginning of the interval A, indicatingthat the increase in gamma-rays detected during the neutron burst 54 maybe dominated by epithermal neutron capture gamma-rays 68. If the neutronburst 54 is short enough, and the separation between neutron bursts 54is significantly longer than a characteristic time τ, the count rateduring the interval A can represented by a single exponential buildup.For example, the following relationship may represent the exponentialbuildup of epithermal neutron capture events:

$\begin{matrix}{{I_{tot} = {I_{inel} + {I_{capture} \cdot ( {1 - e^{- \frac{t}{\tau}}} )}}},} & (2)\end{matrix}$

where τ represents a characteristic time related to the epithermalslowing-down time, and τ is the time from the beginning of the burst 54.Based on the relationship determined in step 134, the data processingcircuitry 44 may determine a count of inelastic gamma-rays 58 in step136. In one embodiment, the data processing circuitry 44 may subtractappropriate fractions of the epithermal count rate and the thermal countrate from the gross count rate in the interval A, leaving a correctedcount of substantially only inelastic gamma-rays 58.

Another embodiment of a method for reducing or eliminating an epithermalneutron capture background from a raw gamma-ray count is the following.From a fit to regions B and C of FIG. 8 using exponential or otherfitting functions, the background count rate immediately prior to theburst and immediate after the burst may be determined. A linearbackground defined from these two points can be used to estimate thebackground underlying the inelastic gamma-ray peak. If the epithermalcapture component is small, then this estimate will be approximate tothe total capture gamma-ray background. However, if the epithermalcapture component is larger, then this estimate will more closelyapproximate the thermal capture component. However, the epithermalcapture component is related to the total capture gamma-ray count rateover the entire cycle, including burst-on and burst-off periods.

Each of the above-described manners of correcting for neutron captureevents may benefit from further techniques for epithermal neutroncapture background reduction. As shown in the curve 120 of FIG. 9 at theinterval A, a gross count rate of gamma-rays increases by almost afactor of two from the start of the interval to its finish. This buildupis due largely to a buildup of epithermal neutron capture gamma-rays 68.If the neutron burst 54 were much shorter, the buildup of epithermalneutron capture gamma-rays 68 would be reduced. As such, a greaterfraction of the gross count of gamma-rays would be due to inelasticgamma-rays 58. In an alternate embodiment, the epithermal capturebackground could be reduced if the neutron burst (interval A in FIG. 7)was sufficiently short that the most of the fast neutrons had not hadsufficient time to slow to epithermal and thermal energies and the gammaray signal would be predominantly comprised of inelastic gamma rays. Theburst would be repeated after allowing sufficient time most of theepithermal neutrons to decay. The remaining thermal neutron capturebackground could then be accurately corrected with conventional methods.

While maintaining approximately the same duty cycle, or fraction ofon-time to off-time, the pulse width of the neutron burst 54 may bereduced. In this way, the epithermal neutron capture gamma-raybackground that occurs during the neutron burst 54 may be minimized. Byway of example, the interval A may be reduced to approximately 4microseconds, the interval B be reduced to approximately 4 to 8microseconds, and the period C may be reduced to approximately 8 to 20microseconds. Other combinations of time intervals are possible forreducing the epithermal neutron capture background with a reduced pulsewidth. As such, the combination of time intervals should not be limitedto the examples stated above.

FIG. 13 depicts a plot 138 illustrating how data 40 obtained from thedownhole tool 12 may be used for determining gas saturation and porosityin the formation 50. The plot 138 may represent a manner of performingthe step 96 of the flowchart 86 of FIG. 5. In the plot 138, an ordinate140 represents a logarithm of a ratio of inelastic gamma-ray counts fromthe far gamma-ray detector 28 to inelastic gamma-ray counts of the neargamma-ray detector 26, represented as log(far/near). Such a ratio ofinelastic gamma-ray counts, free from epithermal and capture gamma-raycontamination, may be representative of the density of the formation 50.A ratio of detector counts is used because certain unknown factors suchas the output of the neutron generator cancel out. Alternatively, if theoutput of the neutron generator is sufficiently stable or is otherwisenormalized, an inelastic gamma-ray count from a single one of thegamma-ray detectors 26 or 28 may be used. In either case, the gamma-raycount may substantially include only inelastic gamma-rays 58 as any ofthe corrective techniques described above may have previously reduced oreliminated epithermal neutron capture background. An abscissa 142represents a logarithm of a neutron count from the neutron detector 22,represented as log(³He). While simulated neutron counts from the neutrondetector 22, log(³He), are employed in the instant example, any signalthat is sensitive to hydrogen index could also be used. Examples ofsignals sensitive to hydrogen index include ratios of neutron countsfrom near and far-spaced neutron detectors, ratios of thermal neutroncapture gamma-rays from near- and far-spaced gamma-ray detector, orsigma. It should be appreciated that that this method of plotting aratio of epithermal-background-free inelastic gamma-rays vs. a hydrogenindex signal is different from the prior art plot of Wilson in U.S. Pat.No. 6,207,953, in which an MeV fast neutron signal is plotted vs. agamma-ray ratio. In the prior art plot, the fast neutron signal is notvery sensitive to hydrogen index, because hydrogen has a relativelysmall influence on transport of MeV neutrons. The fast neutron signal isprimarily sensitive to scattering from the rock matrix atoms, as isdisclosed by Stephenson in U.S. Pat. No. 5,532,482. Additionally, thegamma-ray ratio of the prior art plot contains both inelastic gamma-raysand epithermal capture gamma-rays even though the thermal capturegamma-ray component has been removed. As a result, the prior artgamma-ray ratio is sensitive to both the rock matrix and hydrogen in thefluids and is not independent of the fast neutron signal.

The count rate of epithermal neutrons 62 or thermal neutrons 64 detectedby the neutron detector 22, or log(³He), may decrease as the hydrogenconcentration, or hydrogen index, in the formation 50 increases. Asnoted above, the hydrogen index will decrease if either theliquid-filled porosity (“Φ”) of the formation 50 decreases or if thepore fluid density of the formation 50 decreases, as may occur withreplacement of water by gas (gas saturation, or “S_(g)”). Thus, acombination of measurements primarily sensitive to hydrogen index alonecannot distinguish between these two effects. It should be noted thatmeasurements that purport to measure quantities that are independent ofhydrogen index, such as “a far/near inelastic ratio” of raw GR countsfrom the burst-on (A) timing gate (IRAT), may often be so contaminatedby neutron capture events that they have essentially the same characteras a hydrogen index measurement.

To separate porosity (“Φ”) and gas saturation (“S_(g)”), a secondmeasurement that is truly independent of hydrogen index may be employed.In open-hole, a bulk density measurement may be used, which may bereplaced by the log(far/near) measurement plotted in FIG. 13 Bulkdensity may correspond primarily to the density of atoms in the rockmatrix of the formation 50, which may be complementary to a hydrogenindex measurement. Accordingly, bulk density (log(far/near)) and neutronporosity (log(³He)) may be crossplot to determine true porosity Φ andgas saturation S_(g).

To determine porosity Φ and gas saturation S_(g) using the plot 138 ofFIG. 13, the (log(far/near) corrected inelastic GR count is plottedversus the log(³He) count. The ratios are normalized to the value atzero porosity units (p.u.), such that when the logarithm of the ratio istaken, the 0 pu value is at zero, representing a unity point 144 in theplot 138. Data points 146 are plotted for water-filled pores, whichcontain no gas, at varying porosities Φ. Similarly, data points 148 areplotted for empty pores at the varying porosities Φ. Dashed lines 150represent increasing 152 or decreasing 154 values of gas saturationS_(g) between water-filled pore data points 146 and empty pore datapoints 148 at the listed porosities Φ. At any point along a given dashedline 150, gas saturation S_(g) decreases approximately linearly from 1.0at the corresponding empty pore data point 148 to 0.0 at thecorresponding water-filled data point 146. At any point along the dashedline 150, both the porosity Φ and gas saturation S_(g) may be uniquelydefined.

As should be clear, any point within the boundary defined by thewater-filled pore data points 146 and empty-pore data points 148represents a unique combination of porosity Φ and gas saturation S_(g)values. The transform between points in log(ratio) and log(³He) space toporosity Φ and gas saturation S_(g) space can be calibrated usingformations of known porosity Φ and gas saturations S_(g) or can becalculated using a Monte Carlo technique. The transform can includesmoothly varying functions to connect two spaces or can simply implementa lookup table that interpolates between calibration data. In eithercase, the above-described technique may represent a manner of directlyextracting porosity Φ and gas saturation S_(g) from the raw dataobtained by the downhole tool 12.

With reference to FIGS. 4-13, the downhole tool 12 may be employed indetermining porosity Φ and gas saturation S_(g) by obtaining counts ofinelastic gamma-rays 58. Certain other applications, such ascarbon/oxygen (C/O) logging or the measurement of other elements in theformation 50, may involve gamma-ray spectroscopy involving the inelasticgamma-rays 58. For such gamma-ray spectroscopy techniques, the gamma-raydetectors 26 and/or 28 may include scintillation crystals 32 oflanthanum-halide, such as LaCl. FIG. 5 86 describes a method todetermine porosity and gas saturation for a given set of environmentalconditions.

FIG. 14 illustrates a workflow of an alternate embodiment that accountsfor various environmental effects that can affect the measurement of aformation density sensitive measurement, such as epithermal and thermalcapture corrected far/near inelastic GR counts, and the hydrogensensitive measurement, such as log(³He) counts or far/near capture GRcounts. FIG. 14 is similar to the flowchart 86 of FIG. 5, but addsadditional steps that account for varying environmental effectsincluding, but not limited to; borehole size, casing diameter, casingmaterial, casing thickness, casing position in the borehole, cementcomposition, borehole fluid properties, formation fluid properties,formation gas properties, borehole water salinity, formation watersalinity and formation lithology. All of these environmental variablescan affect the measurements 94 and 162. However, these environmentalvariables can be determined either through another measurement 164 ofsome aspect of the measured gamma-ray and/or neutron from which theenvironmental variable can be determined or by assuming a certainenvironment using external knowledge 166. Once all the environmentalvariables are determined, a proper response function 168 for a givenenvironment can be used, which would be similar, but not necessarily thesame as FIG. 13 and the porosity and gas saturation can be determined170.

Some possible pulsed neutron measurements that are described in priorart that could be used to determine some of the environmental variables162 are (but not limited to): gamma-ray and/or neutron decay rate can beused to determine borehole salinity, formation salinity and/or shalevolume; inelastic gamma ray spectroscopy can be used to determineborehole fluid composition, formation fluid composition and/or lithologyin certain conditions; capture gamma-ray spectroscopy can be used todetermine lithology.

Some possible external measurements 166 that are described in prior artthat could be used to determine some of the environmental variables are(but are not limited to): cement bond logs to determine cementcomposition; production logs to determine borehole fluid properties;pipe integrity logs such as ultrasonic logs to determine casing diameterand thickness; natural gamma ray logs to determine shale volume; openhole logs to determine lithology and other petrophysical properties;pressure measurements to determine borehole and formation fluidproperties and gas densities.

While only certain features of the invention have been illustrated anddescribed herein, many modifications and changes will occur to thoseskilled in the art. It is, therefore, to be understood that the appendedclaims are intended to cover all such modifications and changes as fallwithin the true spirit of the invention.

What is claimed is:
 1. A method comprising: emitting neutrons into asubterranean formation from a downhole tool, such that inelasticgamma-rays are produced due to inelastic scattering in the subterraneanformation; detecting a portion of the inelastic gamma-rays that scatterback to the downhole tool to obtain an inelastic gamma-ray signal,wherein the inelastic gamma-ray signal is substantially free of thermaland epithermal neutron capture background; and determining a property ofthe subterranean formation based at least in part on the inelasticgamma-ray signal.
 2. The method of claim 1, wherein detecting theportion of the inelastic gamma-rays comprises detecting a count of theinelastic gamma-rays to obtain the inelastic gamma-ray signal.
 3. Themethod of claim 1, wherein the portion of the inelastic gamma-rays aredetected using a scintillation crystal consisting essentially of one ormore elements with a resonance integral of less than or equal toapproximately
 12. 4. The method of claim 3, wherein the portion of theinelastic gamma-rays are detected using a scintillation crystalcomprising YAP, BGO, or a lanthanum-halide, or any combination thereof.5. The method of claim 1, wherein determining the property of thesubterranean formation comprises determining, in a processor, a value ofporosity of the subterranean formation, a value of gas saturation of thesubterranean formation, or values of both porosity and gas saturation ofthe subterranean formation.
 6. A downhole tool comprising: a neutronsource capable of emitting neutrons into a subterranean formation suchthat the neutrons inelastically scatter to produce inelastic gamma-rays;and a scintillation crystal capable of detecting the inelasticgamma-rays; wherein substantially all components of the downhole toolwithin a predetermined distance from the scintillation crystal consistessentially of one or more elements with a resonance integral of lessthan approximately
 4. 7. The downhole tool of claim 6, wherein thescintillation crystal consists essentially of one or more elements witha resonance integral of less than or equal to approximately
 12. 8. Thedownhole tool of claim 7, wherein the scintillation crystal comprisesYAP.
 9. The downhole tool of claim 7, wherein the scintillation crystalcomprises a lanthanum halide material.
 10. The downhole tool of claim 6,wherein the scintillation crystal is contained in a scintillator housingconsisting essentially of one or more elements having a resonanceintegral of less than or equal to approximately
 1. 11. The downhole toolof claim 6, wherein the scintillation crystal is contained within athermal neutron shield.
 12. The downhole tool of claim 11, wherein thethermal neutron shield comprises ⁶Li.
 13. The downhole tool of claim 6,comprising a pressure housing that contains the neutron source and thescintillator, wherein the pressure housing consists essentially of oneor more elements having a resonance integral of less than or equal toapproximately
 3. 14. The downhole tool of claim 6, comprising a neutrondetector configured to detect epithermal neutrons.
 15. The downhole toolof claim 6, comprising a ³He neutron detector having a thermal neutronshield.
 16. The downhole tool of claim 15, wherein the thermal neutronshield contains cadmium.
 17. The downhole tool of claim 16, wherein thethermal neutron shield comprises gadolinium.
 18. A system comprising: adownhole tool comprising: an electronic neutron source configured toemit a periodic burst of neutrons into a subterranean formation; and ascintillation detector configured to detect a count of gamma-rays overtime while and after the electronic neutron source emits the periodicburst of neutrons; and data processing circuitry configured to receivethe count and to determine a subset of the count that includessubstantially only inelastic gamma-rays produced by inelastic scatteringof the emitted neutrons with the subterranean formation.
 19. The systemof claim 18, wherein the data processing circuitry is configured todetermine the subset of the count based at least in part on anexponential function fitted to at least two time intervals of the countcorresponding to times when the electronic neutron source is notemitting the periodic burst of neutrons.
 20. The system of claim 18,wherein the data processing circuitry is configured to determine thesubset of the count based at least in part on an exponential functionfitted to a time interval of the count corresponding to a time when theelectronic neutron source is emitting the periodic burst of neutrons.21. A method comprising: receiving, in a processor, a first signalproportionate to a detected quantity of inelastic gamma-rays produced byinelastic scattering of neutrons in a subterranean formation; receiving,in the processor, a second signal proportionate to a hydrogen index ofthe subterranean formation; and determining, in the processor, a valueof porosity of the subterranean formation, a value of gas saturation ofthe subterranean formation, or values of both porosity and gassaturation of the subterranean formation, based at least in part on thefirst signal and the second signal.
 22. The method of claim 21, whereinthe first signal comprises an inelastic gamma-ray count havingsubstantially no epithermal neutron capture background.
 23. The methodof claim 21, wherein the first signal comprises a ratio of inelasticgamma-ray counts detected by a far scintillation detector and a nearscintillation detector, wherein the near scintillation detector isdisposed more closely to a source of the neutrons than the farscintillation detector.
 24. The method of claim 21, wherein the firstsignal comprises a logarithm of count rates of inelastic gamma-rays. 25.The method of claim 21, wherein the second signal comprises a count rateof neutrons detected from a neutron detector.
 26. The method of claim21, wherein the second signal comprises a signal proportionate to acount of a detected quantity of thermal neutron capture gamma-rays,epithermal neutron capture gamma-rays, or both thermal neutron capturegamma-rays and epithermal neutron capture gamma-rays.
 27. The method ofclaim 21, wherein the second signal comprises a signal proportionate tothe ratio of a count of a detected quantity of thermal neutron capturegamma-rays, epithermal neutron capture gamma-rays, or both thermalneutron capture gamma-rays and epithermal neutron capture gamma-raysfrom a first gamma-ray detector and from a second gamma ray detector,wherein the two detectors are at two different axial spacings from theneutron source.
 27. A method comprising: emitting a periodic burst ofneutrons into a subterranean formation using an electronic neutronsource; detecting a raw count of gamma-rays using one or morescintillation detectors, wherein the raw count of gamma-rays includesinelastic gamma-rays caused by inelastic scattering of the neutrons inthe formation and epithermal neutron capture gamma-rays caused byepithermal neutron capture; detecting a count of epithermal neutrons apredetermined distance from the one or more scintillation detectorsusing a neutron detector; and reducing the raw count of gamma-rays by anumber of gamma-rays proportionate to the count of epithermal neutronsto obtain an inelastic gamma-ray count.
 28. The method of claim 27,wherein the count of epithermal neutrons is detected using the neutrondetector, wherein the neutron detector comprises a ³He detector shieldedfrom thermal neutrons.
 29. The method of claim 27, wherein the count ofepithermal neutrons is detected using the neutron detector, wherein theneutron detector is a ³He detector not shielded from thermal neutrons.30. The method of claim 29, wherein detecting the count of epithermalneutrons comprises detecting a count of thermal neutrons and epithermalneutrons using the neutron detector when the periodic burst of neutronsis being emitted, detecting a count of substantially only thermalneutrons using the neutron detector when the periodic burst of neutronsis not being emitted, and reducing the count of both thermal andepithermal neutrons by the count of substantially only thermal neutrons